Giant amplification of interfacially driven transport by hydrodynamic slip: 
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We demonstrate that "moderate" departures from the no-slip hydrodynamic boundary condition 
(hydrodynamic slip lengths in the nanometer range) can result in a very large enhancement - up to 
two orders of magnitude- of most interfacially driven transport phenomena. We study analytically 
and numerically the case of neutral solute diffusio-osmosis in a slab geometry to account for non- 
trivial couplings between interfacial structure and hydrodynamic slip. Possible outcomes are fast 
transport of particles in externally applied or self-generated gradient, and flow enhancement in nano- 
or micro-fluidic geometries. 

PACS numbers: 68.08.-p,68.15.+c, 47.61.-k, 47.63. Gd 
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Introduction - The advent of " microfluidics" and 
" nanofluidics" has motivated the current great interest 
in understanding, modeling and generating motion of liq- 
uids in artificial or natural networks of ever more tiny 
channels or pores Q . Because of the huge increase in hy- 
drodynamic resistance that comes with downsizing, two 
avenues for moving efficiently fluid at such scales have 
been revisited. Both rely on phenomena originating at 
the solid/liquid interface to take advantage of the in- 
crease of surface to volume ratio. 

The first one is the generation of flow within the in- 
terfacial structure by application of a macroscopic gradi- 
ent. Electro-osmosis, i.e. flow-generation by an electric 
field, is the best known example which is commonly used 
in microfluidics 0- But other surface-driven phenomena 
fall in the same category, such as diffusio-osmosis and 
thcrmo-osmosis where gradients of solute concentration 
and of temperature are used to induce solvent flow 0, 0] . 
Their phenomenology is usually best described by an " ef- 
fective slip" velocity, which quantifies the motion of the 
fluid with respect to the solid due to shearing forces in 
the usually thin interfacial layer Q . 

The second is the amplification of pressure-driven flow 
by surfaces such that the fluid hydrodynamically "slips" 
on the solid, as usually quantified by the so-called slip 
length b Q (the distance within the solid at which the 
flow profile extrapolates to zero). Recent efforts in this 
domain have concluded that with a clean " solvophobic" 
surface chemistry one can reach slip lengths up to a 
few ten nanometers 0, but not much more unless to- 
pographic structures are specifically engineered Q. 

In this paperjgeneralizing a point recently made for 
electro-osmosis 0, H, @ , we argue that these two strate- 
gies can actually be synergetically combined, yielding 
strongly enhanced interfacial driven flows on "solvopho- 
bic" surfaces. More quantitatively, we argue that an ac- 
tual " hydrodynamic slip" increases the " effective slip ve- 



locity", which controls all manifestations of the interfa- 
cially driven-phenomena, by a factor (1 +b/L), where b 
is the hydrodynamic slip length and L a measure of the 
interfacial thickness. This ratio can be of order ten to 
hundreds in realistic situations, so that the enhancement 
described here can be very large. This synergy may lead 
to more efficient transduction of electrical, chemical or 
thermal energy into mechanical work in micro-devices. 

Beyond the nano/microfluidic interest in moving flu- 
ids in tiny solid structures, our considerations also ap- 
ply to the reciprocal interfacially driven motion of solid 
particles in solution. We thus predict enhancement of 
electrophoresis, diffusiophorcsis, and thcrmophoresis (in- 
duced respectively by gradients in electric potential, con- 
centration of solutes and gradients of temperature) when 
solvophobic particles are dispersed in solution. Our anal- 
ysis may also be of relevance to the "swimming" of ar- 
tificial or natural organisms by self-generation of such 
gradients 0, E3, EH H3 

To exhibit the physics at work, we first focus on 
diffusio-osmosis with a single neutral solute species, 
in the simplest geometry of a flat uniform interface. 
Using a continuum description for hydrodynamics 
with slip, we derive the (1 + b/L) enhancement factor 
for that situation. A formal generalization to other 
intcrfacially-driven phenomena is then presented. Fur- 
ther, a molecular dynamics study of diffusio-osmosis in 
a thin slab geometry quantitatively conforts the picture. 
We end with a brief discussion of the case of charged 
solutes (in particular electro-osmosis) and of the motion 
of finite-size particles. 

Consider a flat homogeneous solid surface y(x, z) = 
0, with an incompressible liquid of bulk viscosity r/o 
in the y > space. Slip is decsribed through the 
Navier boundary condition (BC) for the velocity field 



bd, 



yVi\y= 



i>i\ y= o for i = x, z, with b the distance 



in the solid at which the linearly extrapolated velocity 
becomes zero (see Figure 1). In a slightly more general 
approach the hydrodynamic "weakness" of the interface 
shows up in a y-dependent viscosity rj(y), while requiring 
v\ y= o = 0. The Navier BC is recovered using the ansatz 
V~ 1 (y) = Vo 1 ^ + b.6(y)), sketching slip in terms of a 
very thin vacuum layer of very low viscosity "between" 
the liquid and the solid. 

Diffusio- osmosis for a neutral solute - Suppose that 
the solution contains a single neutral solute, at bulk 
concentration cq, which interacts with the wall through 
a short-range potential U(y). In the dilute limit, at 
equilibrium the distribution of the solute is c eq (y) = 
Co exp(— fc^p )■ If a concentration gradient dco/dx is 
applied along x over long distances (compared to the 
range of the potential), equilibration of concentration 
and pressure is fast along y (compared to the relaxation 
time of the gradient), so that c ~ Cq(x) exp(— and 
—d y p(x,y) + c(x,y)(—d y U) = 0. This leads to the "os- 
motic" equilibrium p{x, y) — 1cbTc(x, y) = po — fcsTco(x) 
|l3j . with po the constant bulk pressure. As a conse- 
quence a pressure gradient along x sets in (only) within 
the thin interfacial layer, d x p = kBTd x (c — Co), which 
generates shear there through the hydrodynamic balance: 
—d x p{x, y) + d y (rj d v v x ) = 0. The fluid velocity increases 
accordingly through the interfacial layer to reach a finite 
value v s , the "effective slip velocity" of the liquid past 
the surface due to the applied gradient along x (Figure 1 
sketches the case of a solute attracted to the wall, T > 0). 
Integrating twice along y and using the Navier BC: 

v s = -(k B T/ m )TL(l + b/L).^ (1) 

where T = J °°dy[e~ u ^ v ^ kBT — 1] is a length measur- 
ing the excess of solute in the vicinity of the surface 
(U is positive and T negative for depletion), and L = 
T" 1 f °°dyy{e~ u ( y ^ kBT — 1] measures the range of inter- 
action of the potential. Equation Q is the classical for- 
mula of Anderson and Prieve |l3j , times the amplification 
factor 1 + b/L : this quantifies how hydrodynamic slip al- 
lows to generate a larger " effective slip" v s away from the 
surface (Figure 1). Physically v s results from the balance 
between viscous shear stress at the interface corrected for 
slip, rjov s /(L + b), and the (integrated) body force within 
the interface layer : — -^(rksTco). 

The slip induced enhancement can actually be very 
large. For molecular interactions between neutral solutes 
and a solid L is very small, e.g. ~ 0.3 nm, so with 
b ~ 20 — 30 nm for water on hydrophobic substrates 
[HEli the amplification factor can be up to 100 ! 

Formal general argument - We now generalize this re- 
sult. For a generic interfacial structure, denote a n and a t 
the stresses normal and tangential to the interface which 
develop in a thin layer close to the solid Q. At equilib- 
rium, the situation is invariant by translation along x, so 
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FIG. 1: Diffusio-osmosis for a neutral solute attracted 
to the solid (in grey) . The larger bulk concentration on 
the left results in a higher accumulation of solute at the 
interface (thickness L) "squeezing" the fluid against the 
wall. A pressure gradient results in this L-thick interfacial 
layer, which induces shear there and a flow opposite to the 
concentration gradient. In contrast with the no-slip case 
(b — 0, flow depicted by the dot-dashed line), hydrodynamic 
slip (i.e. a non-zero extrapolation slip length b) allows these 
stresses to generate a larger "effective" diffusiophoretic slip 
i) s (to be distinguished from the very local slip velocity right 
at the wall). 

"V» = o~n{y) an d c t = o~ t (y), and the hydrostatic pressure 
p(y) is determined by force balance —d y p + d y a n = 0, 
-d x p + d x a t = 0, yielding p(y) = p + a n (y), with again 
Po the constant bulk pressure. 

If a small far-field gradient of an observable O (concen- 
tration, potential, temperature) is applied along x then 
the interfacial stresses vary slowly along x too. Pressure 
equilibration is fast in the y direction and —d y p + d y o- n = 
yields p{x,y) = po + a n {x,y). The resulting lateral 
imbalance of pressure, within the interfacial layer, gen- 
erates shear along x as described by the force balance 
—d x p + d x a t + dy(n(y)d y v x ) = 0. Again, this gener- 
ates an "effective slip" v s which reads (using r](y) and 
v(y = 0) = 0): 
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with Ti(x,y) = a n — a t the interfacial stress anisotropy. 
If the structure of the interface varies slowly along x, 
E ~ S e9 (y, 0(x)) where 0(x) is the "outer" value of the 
field O outside the interfacial layer, so that the effective 
slip generated by dO/dx is 
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The integral in the bracket quantifies the specific con- 
tribution of the hydrodynamic slip. For a slip length b 
[using f]o/j](y) = 1 + b8(y)], we obtain our main result: 
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where T = J °°dy.^ff-{y), and L = T 1 j^dy.y.^iy) 
is a measure of the thickness of the stress-generating in- 
terface, that depends on the observable O considered. 
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The case of (neutral solute) diffusio-osmosis is recov- 
ered with: a n — fcsT(c — Co), a t = 0, O = ksTco and 

^%ff- = (c eq - c )/c = e ^ - l. 

The results obtained so far rely on a continuum 
description of the interface hydrodynamics. To demon- 
strate that the enhancement persists in a more realistic 
context, we turn to a slab geometry, that we will analyze 
using numerical simulations. 

Diffusioosmosis in a channel - Let us consider a chan- 
nel of width H . In the linear response regime a symmetric 
matrix relates the fluxes (per unit length) in the x direc- 
tion (Q,J) to the gradients (— Vtt, — V/z) that generate 
them, with Q the total flow rate, J the total solute cur- 
rent, 7r = p — ksTc the pressure corrected for osmotic 
effects, and /j, ~ un + A;bT1ii(c) the chemical potential 
of the solute [Til Tlij . Equivalently, diffusio-osmosis is 
best decribed by the following matrix M quantifying net 
transport through the channel: 



(5) 



Onsager reciprocity relations require A/12 = A/21 
which we explicitly checked by solving the continuum 
hydrodynamic problem with a slip length b on the 
two walls in the two situations Vp 7^ 0, Vco = and 
Vp = 0, Vcq ^ 16] . hi the latter situation and for 
channels wider than the interfacial structures (H 3> L), 
we obtain as expected a plug-like flow driven by a con- 
centration gradient : Q = A/i 2 (— /csTVc /c ) = H.v s 
with v s the slip velocity given in Eq. and 
A/12 = ""^(i + k)- For thinner (nano) channels, the 
overlap between the interfacial layers must be taken into 
account |l6j . 

Numerical Simulations - We then conduct Molecu- 
lar Dynamics simulations of a fluid system composed 
of solvent+solute particles, confined between two par- 
allel solid walls composed of individual "solid" parti- 
cles fixed on a fee lattice |l7|. Interactions between 
the three types of particles are of Lennard- Jones type, 

"(7) 12 -^(7) 6 " 
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with identical intcr- 
and molecular diameters a (a, (3 £ 



E7a/j(r) = 4e 
action energy e 
{solute, solvent, walls}). Tuning the parameters u a ^ we 
can vary (i) the wettability of the solvent on the wall by 
tuning u S oivcnt,waii, and (ii) the relative attraction or de- 
pletion of the solute to that wall (by tuning M W aU,solute 
for a fixed u wa ii, solvent)- Periodic boundary conditions 
are used along x and z (box size l x = l z = 16a), and the 
inter-wall distance is l y = H = 20.8a. Temperature is 
kept constant by applying a Hoover thermostat to the z 
degrees of freedom (i.e. perpendicular to flow and con- 
finement). Solvent density is p/cr 3 ~ 0.9, and bulk solute 
concentration coa 3 ~ 0.02. Rather long runs (~ 5.10 6 
timcsteps) are performed to obtain good statistics. 
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FIG. 2: Cross coefficient M12 for a slab geometry (normalized 
by a no-slip reference M^, = —oHY), against normalized in- 
terfacial enrichment in solute T/a, for wetting to non-wetting 

solvents W W all, B olvent = 1-0 (a), 0.5 ("),0.3 (•), and Uwall.solute 

in the range [0.1 — 1.1]. The dashed lines are the theoretical 
prediction M 12 /M$ = (L + b)/a ~ b/a using the measured 
slip length b which is significant for the non-wetting cases: 
b ~ 12 — 16a (■) and b ~ 20 — 10a (•). For large positive ad- 
sorption, the enhancement decreases as does the slip length, 
because the adsorbed solute increases the fluid/solid wetting. 

To determine the cross coefficient A/12 = A/21, the 
most efficient route is to apply an external volume force, 
fo = — Vp, to the fluid in the x direction, so as to gen- 
erate a pressure-driven flow. We then measure the so- 
lute excess current, J — cqQ, associated with the con- 
vective motion of the solute 01 1 an d obtain M 2 i = 
(J — CqQ)/(— Vp), according to eq. JSJ (we check lin- 
earity of the reponse to the external force) . Eventually 
we extract the adsorption length T from the equilibrium 
solute density profile c(y) as T = | J slab dy[c(y) / cq — 1]. 
To narrow our exploration, we focus on the ideal solution 
of solvent and solute molecules identical but for their 
interactions with the walls. We take "solvent, solvent = 

"solute, solute = "solvent, solute = 1-2, and Consider tlu'CC 

solvent-wall situations "wall, solvent = 0.3,0.5,1.0 (going 
from non- wetting to wetting) and solute- wall interactions 
"wall, solute hi the range [0.1,1.1]. In all cases, the hy- 
drodynamic velocity profiles are parabolic, which allows 
us to extract the viscosity r\ and the slip length b [l9j . 
In agreement with previous work [l9j and experimen- 
tal results, slip is significant for a non-wetting solvent 
(6 - 20 - 40a for it wa ii, solvent ~ 0.3 - 0.5), and negligible 
for a wetting solvent (b < a for u wa ii, solvent = 1)- 

Fig. |21 displays the outcome of our simulations for 
the cross coefficient A/21 = M12, normalized by a value 
M^2 = ^-aHT, which corresponds to a reference sit- 
uation with a no-slip BC and fixed L = a. In line 
with our theoretical arguments, Mi 2 is strongly enhanced 
- M 12 /M 1 ( 2 ) up to 40 here- for non-wetting solvents 
("wail, solvent = 0.3,0.5, top data), i.e. systems with slip 
lengths of a few tens of molecular diameters (i.e. roughly 
~ 5 — lOnm). On the other hand, Mi 2 is much smaller 
(A/12 ~ M}^) for a wetting solvent (" wa ii, solvent = 1, 
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bottom data), associated with a no-slip BC. 

More quantitatively our MD results compare suc- 
cesfully to the theoretical prediction rewritten as 
M 12 /M[f = (L + b)/a ~ b/a (since L < a < b), see 
Figure [21 provided we use the slip length b extracted 
from the simulation for each case. In particular, the am- 
plification decreases for large positive adsorptions T > 0, 
due to the decrease of the slip length : accumulation of 
"wetting" solute at the solid- liquid interface reduces the 
effective "solvophobicity" [2(|- For a depleted solute, 
(r < 0), this "saturation" effect is essentially absent (b 
is nearly constant) allowing for large enhancements. 

Electrolyte solutions: electro- and diffusio- osmosis - 
For charged solutes, the interfacial structure is the elec- 
trical double layer, of typical thickness the Debye length 
0, usually in the nanometer range (1 — 30 nm), 
and we anticipate L ~ k . The enhancement of in- 
terfacially driven phenomena (electro-osmosis, diffusio- 
osmosis, thermophoresis) over solvophobic surfaces {b in 
the 20 — 30 nm range) should thus be somewhat smaller 
than for the neutral solute case, but still significant. 



As a check of L 



we incorporate a finite slip 



length b in the usual description of these phenomena 
and compute the enhancement factor (1 + b/L) in 



j/d»l»= ' 



equation (4) For electro-osmosis L = 

with <j) eq (y) the equilibrium electrostatic potential in the 
double layer, so for weakly charged surfaces L ~ 
in agreement with 0, Ej- For diffusio-osmosis and a 
1:1 electrolyte, wc obtain a more complex formula, with 
L ~ k _1 /2 for weakly charged surfaces. 

Transport of particles - All the above applies to the re- 
ciprocal motion of particles in concentration or potential 
fields, in a way that can be quite directly quantified pro- 
vided the surface is locally flat and homogeneous at the L 
scale following 0, 0, El • Classically for interfacial driven 
effects, considering finite-size objects such as a spherical 
particle of radius n>I, allows one to discuss the pos- 
sible feed-back of the generated flow on the interfacial 
structure where it originates |l3j . We compute here the 
diffusiophorcsis of such a sphere generated by a steady 
background gradient of neutral solute, adapting the clas- 
sical no-slip analysis of [P3T | . Including hydrodynamic slip 
(non-zero 6) enhances the surface/ liquid effective slip as 
described by (1), but also the convection of solute in the 
interfacial region, which affects the steady-state concen- 
tration field of the solute (diffusion coefficient D) around 
the particle. We find the velocity U of the particle in a 
solute gradient Vco to be: 



U 



k B T 



.LT.- 



l + b/L 



l + (l + {v + b/L)Pe)(T/a) 



•Vc (6) 



shape of the potential. For moderate values of Y <C a, 
the usual slip enhancement factor (1 + b/L) prevails, 
and for b ^> L the formula reads U ~ (k B T / 'r)o) .bY .V c . 
For Vc - 10- 3 mol/cm 4 , Y ~ L ~ A, and b ~ 30nm, 
this leads to U > /Um/s (in contrast to ~ 5nm/s for the 
no slip case !), comparable to experimental observations 
of chemical self-propulsion F° r smaller particles or 
stronger solute adsorption (Y/a » 1), the effect of slip 
saturates for large b/L (the large generated flow "erases" 
partly the original interfacial gradients), with a max- 
imal velocity U max = (fcsT/ryo).p^.Vco independent of b. 



Conclusion - Hydrodynamic slip can very significantly 
enhance many interfacially-drivcn phenomena on smooth 
"solvophobic" surfaces. This is of relevance for the trans- 
port of fluids in small channels, and of particles in solu- 
tions. A related target is the modelling and engineering 
of the self-transport of chemically-driven swimmers, for 
which the hydrophobicity of the surface is thought to 
play an important role |l2j. Further study is necessary 
to go beyond the model smooth surfaces considered here, 
so as to assess the effect of topographic or chemical het- 
erogeneities at various scales (e.g. roughness can poten- 
tially cither increase or decrease slip effects in channels 
depending on whether or not it leads to gas entrapment). 



with Pe = (fcfjT / Dr^o)LTcg and v a dimcnsionless quan- 
tity of order 1 defined in |l3| that depends on the exact 
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